Astronomy & Astrophysics manuscript no. vida'et'al 


© ESO 2009 


September 28, 2009 





o 
o 

(N 

oo 

00 
(N 

P? 

00 

43 
6 

o3 



> 



o 
o 



13 



Photospheric and chromospheric activity in V405 And 

An M dwarf binary with components on the two sides of the full convection limit 

K. Vida ' 2 , K. Olah 2 , Zs. Kovari 2 , H. Korhonen 3 , J. Bartus 4 , Zs. Hurta 1,2 , and K. Posztobanyi 5 

1 Eotvos Lorand University, Department of Astronomy, H-1518 Budapest, PO Box 32, Hungary 
e-mail: vidakris@elte.hu 

2 Konkoly Observatory of the Hungarian Academy of Sciences, 1525 Budapest, PO Box 67, Hungary 

3 European Southern Observatory, Karl-Schwarzschild-StraBe 2, 85748 Garching bei Miinchen, Germany 

4 Astrophysical Institute Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany 

5 AEKI, KFKI Atomic Energy Research Institute, Thermohydraulic Department, H1525 Budapest 1 14, PO Box 49, Hungary 



Preprint online version: September 28, 2009 



ABSTRACT 



We investigate the fast rotating (P orh = P rol = 0.465a! ) active dwarf binary V405 And (M0V+M5V) using photometric BV(RI) C 
and optical spectroscopic data. The light variation is caused by the combined effect of spottedness and binarity with a small eclipse. 
From the available light and radial velocity curves we estimate the system parameters. Three flare events happened during the obser- 
vations: two were found in the spectroscopic data and one was observed photometrically in BV(RI)c colours. An interesting eruptive 
phenomenon emerged from the photometric measurements which can be interpreted as a series of post-flare eruptions lasting for at 
least 3 orbits (rotations) of the system, originating from trans-equatorial magnetic loops, which connect the active regions in the two 
hemispheres. The two components of V405 And have masses well over and below the theoretical limit of full convection. This rare 
property makes the binary an ideal target for observing and testing models for stellar dynamo action. 
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1. Introduction 

Low mass stars are in the focus of studying stellar structure: 
these kind of objects make up the majority of the Galaxy, yet, 
because of their low luminosity, it is not easy to study them in 
detail. Low mass stars with fast rotation show magnetic phe- 
nomena (spots, plages, flares, strong X-ray emission) due to 
some kind of magnetic dynamo in their interiors. Stars with 
radiative core and convection zone are thought to operate aQ. 
dynamos. The a-effect creates poloidal field from toroidal one 
by helical turbulence in the convection zone, whereas £2-effect 
produce toroidal field from poloidal one by differential rota- 
tion. The place of the dynamo process is thought to be a thin 
layer between the core and convection zone, called tachocline, 
where the shear is the strongest. Below about ~O.32M stars are 
thought to be fully convective, yet show all known magnetic phe- 
nomena as well. The source of the magnetic field in these low 
mass stars could be a 2 dynamo generating strong, long-lasting, 
axisymmetric magnetic fields by turbulence, or by some other 
mechanism involving differential rotation. However, the mass 
limit (or the spectral limit of M3.5) of the full convection is 
not well established, since strong magnetic fields may shift this 
boundary towards lower masses ( Mul lan & MacDonaldl [200 1 , 
Chabrier et al. 2007). For low mass (M < 0.8) active stars a sig- 
nificant discrepancy is found between models and observations 
of stellar radii, i.e., theory predicts ~10% lower rad ii, while tem - 
peratures are underestimated by about 5% (see e.g. Ribas 2006). 

An interesting approach is to study low mass stars above and 
below this hypothetical mass limit together in close binary sys- 
tems. As to our knowledge, only two such systems are known. 



2MASS J04463285+ 1901432, in NGC 1647 dHebb et al.| [2006) 
has an orbital period of 0.619<i and masses of 0.47 and 0.1 9M , 
however, its V - 19.34 prevents more detailed investigation of 
this binary. 

The binary studied in this work is an X-ray bright ob- 
j ect, listed in the R OSAT All-Sky-Survey as RX J0222.4+4729 
dVoges et alj 119961) . later named as V405 And classified as 
RS CVn type in the G eneral Catalogue o f Variable Stars 
dKazarovets et alj 120001: ISamus et alj I2009I). The on l y de- 
tailed study of V405 And is from lChevalier & Ilovaiskvl d 19971) 
(CI97 hereafter), who find that the binary has an orbital pe- 
riod of P = QA65d showing a small (near grazing) eclipse, and 
present photometric observations in B and V, radial velocity data 
and Ha study. The components have spectral types MOV and 
M5V derived by CI97, both of them showing strong Ha emis- 
sion lines, of which the primary is modulated by the rotation. 
The original classification of the system is BY Dra binary by the 
above authors. Note, that the definitions of the 'RS' and 'BY' 
types overlap in the GCVS, and for V405 And both classifica- 
tions seem to be valid. In CI97 orbital elements and estimates of 
the masses around 0.2 and 0.5M Q , and radii of 0.78 and 0.25- 
0.30R Q are given as well, which were used as starting point for 
our analysis. 



2. Observations and data reduction 

Photometric observations were carried out using the lm RCC 
telescope of Konkoly Observatory at Piszkesteto Mountain 
Station equipped with a Princeton Instruments 1 300 x 1 300 CCD 
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Fig. 1. Finding chart for V405 And from a V-band frame 
of the lm RCC telescope . Numbering is the same as in 
IChevalier & Ilovaiskyld 19971) . 



and the 60cm telescope of Konkoly Observatory at Svabhegy, 
Budapest with a Wright Instruments 750 x 1 100 CCD camera. 

Observations were collected during 28 nights: 13 nights 
with the 60cm telescope (February 8 - October 1, 2007; 
2454140-2454375 HJDs), and 18 nights with the RCC telescope 
(September 28, 2007 - November 18, 2008; 2454372-2454789 
HJDs), on three nights both telescopes were used simultane- 
ously. 

After correcting for atmospheric extinction and transforming 
the data to standard photometric system, we matched the two 
data series by fine-tuning the telescope constants within their 
errors. 

A finding chart is plotted on Fig.Q] V405 And is marked as 
Star 1 (using the numbering of CI97). Two of the comparison 
stars (Star 3, Star 6) used by CI97 turned out to be variable. 
Star 3 shows changes on a longer timescale of 20-30 days. Star 
6, showing a period of ~ 0.658ii, is possibly a W UMa type star. 
This effect would only cause small shifts, if any in CI97, since 
they used the average value of the five comparison stars. 

In this work GSC 03298-00148 (Star 7) was used as a com- 
parison star. Magnitudes of Star 7 were calculated using stars 
2,3,4 and 5 separately and the results were compared. The dif- 
ferences were 0.03 magnitudes at maximum. According to the 
results, we chose finally Star 2 for defining the magnitudes of 
the comparison star. 

As mentioned in CI97, there is a star near V405And, at 
about 4.5 arc-seconds to the south-east, which is blended with 
the variable on most of the frames. We determined the magni- 
tudes of the star near V405 And using frames taken on 2007 
September 29, October 01 and October 02. IRAFQ PHOT rou- 
tine of DAOPHOT package with small aperture, and NSTAR 
PSF-fitting routine were used. There was a systematic differ- 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Table 1. Average magnitudes of stars on V405 And field. 



Star Nn * 


y 


B - V 


V - Rc 


V - Ir- 


1 (V405 And) 


11.25 


1.36 


0.94 


1.86 




11.10" 


1.44 


0.89 


1.74 


2 


12.80 


0.43 


0.28 


0.55 




12.80 


0.43 


0.28 


0.55 


3 


13.60 


0.43 


0.28 


0.57 




13.60 


0.41 


0.28 


0.54 


4 


11.05 


0.44 


0.28 


0.54 




11.05 


0.43 


0.28 


0.54 


5 


12.99 


0.89 








12.96 


0.94 






6 


11.67 


0.68 


0.45 


0.83 




11.65 


0.72 


0.42 


0.79 


7 


10.67 


0.76 


0.47 


0.87 



* The star numbers are the same as in CI97. Star 7 is not in the field of 
view in CI97. 

** Values in italics are from CI97. 

ence (order of 0.01 magnitudes) between aperture and PSF pho- 
tometry (aperture photometry gave brighter magnitudes). Since 
the two stars are separated only by 12.6 pixels, V405 And ex- 
tends into the aperture of the fainter star, so the result of the 
PSF photometry was accepted. The magnitudes of the faint star 
are 16.52 + 0.06, 14.92 ± 0.05 (15.1 in CI97), 13.75 + 0.03 and 
12.13 + 0.03 in B, V,Rc and Iq passbands, respectively. To test 
the correctness of these values, we compared PSF-magnitudes 
of V405 And with magnitudes from aperture photometry. We 
set an aperture large enough to contain both V405 And and the 
faint star, then subtracted the intensity of the faint star from the 
result. These magnitudes from aperture photometry and PSF- 
photometry match reassuringly within a few hundredths of mag- 
nitude. 

Differential aperture photometry of V405 And was done us- 
ing DAOPHOT routine. The aperture was chosen to contain the 
faint star on all frames, since the two stars were not resolved on 
most frames. The intensities of this star were subtracted from 
the light curve after transforming to standard Johnson-Cousins 
photometric system. The resulting light curve is plotted in Fig. [2] 

Interstellar extinction was checked for V405And using 
NASA/IPAC Extragalactic Database (NED), giving a value of 
E(B - V) = 0.182 mag for total extinction. According to CI97, 
the distance of V405 And is 25-35 pc, so the effect of interstel- 
lar extinction would be only a few thousandth magnitude. This 
value is much smaller than the errors of the photometry (~0.01- 
0.02mag), so it can be neglected. 

Altogether 97 high-resolution spectra was retrieved from 
the ELODIE archive a t Observatoire de Haute-Provence (OHP) 
dMoultaka et al.l 120041) . In 1994 38 spectr a were observed: 27 
between October 18-23, and 11 between December 12-16. 
Another 59 unpublished spectra were obtained between October 
25-30, 1998. The spectra from the archive cover the wavelength 
range of 4000-6800A with a resolution (A/AA) of 42000. 

3. Analysis 

Results in Table Q] show that while the average magnitudes of 
our measurements agree quite well with those of CI97 in case 
of Stars 2 to 6, values for V405 And seem to differ significantly. 
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velocity curve from CI97. The binary model light curve is re- 
moved from the observed one using the following equation: 
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Fig. 2. a) Light curve in V passband. The scale is smaller than the 
actual changes of the star caused by a flare at the beginning of 
the observation. bl,b2) Four-colour light curves of V405 And, 
and the fitted light curve solution containing both binarity and 
spottedness. The phased light curves contain observations from 
Julian date intervals shown on the plots. 



V405 And was 0. 15 magnitude fainter during our observations in 
V band than in 1995. This, and the decrease of V—Rq and V—Ic 
indices indicate higher spottedness. In the same time, the bluer 
B — V show the presence of more faculae. Thus, we conclude that 
the activity is enhanced compared to the 1995 level. 



3. 1 . Spots: photometric modelling 

Fo r finding light curve solut ion we adopted the method described 
bylStrassmeier ^tall (|2008|) . The method uses an iterative mod- 
elling for the binarity and spottednes s separately. The b inary 
model was computed using PHOEBE (Prsa & Zwitter 2005) si- 
multaneously for our four-colour measurements and the radial 



pot 



obs 



model 



+ F 



unsp i 



(i) 



where F spot is the flux of the primary changing due to spotted- 
ness, F bs is the observed light curve, F mo d e i is the binary model 
constructed with PHOEBE, and F unsp is the unspotted flux of the 
system (for model parameters see Table [2). F spot is a light curve 
which contains only the changes caused by the spottedness, and 
is free from the effects of binarity (eclipse and changes by the 
distorted shape of the stars ). The true spotted light curve was 
modelled with SpotModeL dRibarik et al. 2003), which uses the 
analytic approach with circular spots of Budding ( 19771) . The cir- 
cumpolar part of the spotted region accounts for the long-term, 
overall brightness change, whereas spotted parts outside this re- 
gion cause the rotational modulation. The latitude of the spots is, 
at least partly, determined by the ratio of buoyancy and Coriolis- 
force. On fast rotating stars the Cor iolis-force is dominant which 
moves the spots to high latitudes dChoudhuri & Gilmanlll987l) . 
While spot models can determine spot longitudes quite reliably, 
photometric observati ons have very low info rmation content on 
spot latitudes (see e.g. lKovari & B artus 1997). Assuming circu- 
lar spots on the surface is an approximation with small number 
of parameters. Sunspots and active nests on the Sun give ex- 
amples of circular active regions. Note, that as 'spots' we mean 
regions with dark spots and bright faculae, the (average) tem- 
perature of the region is made up from the cool and hot parts, 
relative to the photosphere. 

We assumed two spotted regions on V405 And with average 
F; pot -3300K on th e high- latitude regions of the primary. Tests 
of iKovari & Bartusl (1 19971) showed that analytic models using 
two circular spots describe the observed spotted light curves sat- 
isfactorily within the precision of the photometric observations. 
Spot temperature was chosen by modelling a V - Iq light curve, 
and the result was accepted for all passbands. For finding the 
best fitting spotted light curve, the size and position of the spots 
were searched. During our observations the star was more than 
0. 1 magnitudes fainter due to higher spottedness than the faintest 
observed value in CI97, thus we used the brightest part of the V 
and B—V light curves of CI97 as the (faintest possible) unspot- 
ted B and V magnitudes. Unspotted Rc and Iq magnitudes were 
calculated from the color indices. The resulting spotted curve 
was then removed from the observed one, and these steps were 
repeated until we got a satisfactory fit for the spots and bina- 
rity (Fig. [3} separately, and for the observed light curve itself 
(Fig.E>). 

Using this modelling method it is possible to determine the 
physical parameters of the binary (see Table [2j. For initial v alues 
of surface temperatures we u sed the tables of lFlowen (1 19961) and 
IVand enBerg & Cleml d2003l) . The shape of the light curve is not 
very sensitive to the changes in the temperature, thus the fitting 
of these values is quite uncertain. Our result for the mass ratio is 
within lcr from the value that was spectroscopically determined 
inCI97 (q = 0.38 ± 0.04). 

The spot configuration is relatively stable throughout the ob- 
servations, positions and sizes of the spots change only slightly. 
According to the model, one spot is situated on the northern 
hemisphere, and another one on the southern hemisphere about 
150° from each other in longitude. From the modelling results 
the longitudes of spots (i.e., phases) are accurate within a few 
degrees, but latitude results are just approximations. Using the 
supposed r spot we find that about 25% of the stellar surface is 
covered with active regions. The stability of circular spot mod- 
els are discussed in details in Kovari & Bartusl (Q997). With the 
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Fig. 3. Phased BV(RI)c light curves from 2007 after subtracting 
the spot model from the observations. This dataset was used for 
determining the binary system parameters with PHOEBE. 

Table 2. Binary parameters from the system modelling and spot 
model parameters. 



— P rot ~ 


= 0.46543d * 


a = 


= 2.246fi ** 


q -- 


= 0.42 ± 0.01 


i - 


= 66.5 ± 1° 


Mi = 


= 0.49M o 


M 2 -- 


= 0.21M o 


T\ -- 


= 4050 ± 200 K 


T 2 = 


= 3000 ± 300 K 


Ri 


= 0.78 ± 0.02tf o 


R 2 -- 


= 0.24 ± 0.04tf o 


v sin i = 


= 85 km/s* 


e - 


= 0.0* 



^unsp 


= 12.42 


V - 


= 11.00 


-^unsp 


= 10.07 


Ainsp 


= 9.17 


^spot 


= 3300K 



adopted from CI97 
* asin i adopted from CI97 



a 
x 




Fig. 4. Equivalent widths of the Ha line. Data obtained on 
October 18, 1994 are plotted with filled points, empty points 
show measurements from October 21/23 and December of 1994, 
spectra from 1998 are plotted with crosses. Arrows show times 
of flares in 1994 and 1998 (see Section[373]l. 
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Fig. 5. Light curve observed on September 29-30, 2007 at 
Piszkesteto. The binary and spot model is removed from the ob- 
served light curves. Dashed lines show zero level for each pass- 
band. Data from the 29th and 30th September are plotted with 
empty and full circles, respectively. 



help of this simple spot model we could satisfactorily separate 
the light variation originating from binarity and activity. 

3.2. Hot regions in the chromosphere and photosphere 

The chromospheric activity of V405 And was investigated using 
Ho- line profiles. Equivalent widths (Fig. [4]l were measured by 
fitting Gaussian profiles using SPLOT task in IRAF, after the 
continuum was fitted by polynomials and removed. The equiva- 
lent widths vary between 1-2 A showing a small rotational mod- 
ulation. This variation seems to be fairly regular during the ob- 
servations, except for the time of the flares in October 1994 and 
1998 (see Section [3~3"l ). The observed variation indicates perma- 
nent plage structure in the stellar atmosphere. 



In September 29-30, 2007 an interesting phenomenon was 
observed (see BV(RI)c light curves in Fig- Et. During two con- 
secutive nights, at the beginning of the observations the bi- 
nary showed excess brightness in all bandpasses which gradu- 
ally decreased to the unspotted flux. After about one third rota- 
tion the brightening started again (the system configuration dur- 
ing the event is plotted on Fig. |6j. The excess brightness was 
stronger during the second night (between the two sets of ob- 
servations was one full rotation). The ratios of the peak intensi- 
ties in BV(RI)c filters at the beginning of the observations differ 
from that produced by a usual flare event at or near its peak in- 
tensity. 

Assuming that the brightening is caused by hot regions on 
the star we can make a model of the surface using SpotModeL. 
Results show that two hot areas separated in longitude by ~90- 
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Fig. 6. Schematic plot of the eruption times and visibility on 29- 
30 Sept, 2007. Arrow at the bottom shows the line of sight. The 
figure is plotted with appropriate scale. 



120° describe the observed changes well. Making use of two 
passbands, from B—V and V—Ic indices we get hot spot tempera- 
tures of 5 1 30+270K and 5680+300K, respectively, which means 
these regions are at least 1000K hotter than the surrounding area. 
At maximum of the brightness B — V = 1.31 yields 4300K for 
the temperature dFlowerlll996l) . meaning that the whole visible 
stellar surface seemed hotter by ~150 K, i.e.; about 10-12% of 
the visible surface of the star is covered with hot regions. 

An explanation of the phenomenon could be a long-lasting 
(at least 1 .5 days, or 3 rotations) flare event on V405 And. During 
the first night's observations we see only the decreasing part of 
a giant flare that might have occurred during daytime, before the 
observations started. At the beginning of the second night we 
observe a peak which could then be a post-flare eruption. Two 
smaller eruptions happened later that night, the first was small, 
the second was more powerful and started to decrease still before 
the end of the night. 

A similar long-lasting event has been observed on AU Mic . 
This event was described and modelled in lKatsova et al.l {1999). 
The authors, comparing their results with observed solar phe- 
nomena, supposed long-lasting giant post-flare loops with a size 
of 1-2 stellar radii, and energy balance through continuing re- 
connections. 

Transequator ial magnetic loops are common features on the 
Sun. lHarra et a il (120031) studied the energy mechanism of such 
a solar loop system and found brightenings and flare-like events 
during 2 days of observations. On V405 And from photomet- 
ric observations we suggest 1-1 dominant active region in both 
hemispheres of the primary separated by ~150° in longitude. 
We could explain the observed phenomenon on September 29- 
30, 2007 with a transequatiorial loop system which causes the 
brightenings and the flare-like events. As the star rotates, the 
emitting part of the loop system eventually turns in and out of 
view. This scenario would explain the turn-backs of the light 
curve seen, e.g., in B band at 0.3 phase (Fig. |5), and also the 
intensity ratios in different filters, which are comparable to the 
post-flare event of Fig. [7] 



3.3. Flares 

During the observations, three noteworthy flare events happened. 
Two of them were observed with the ELODIE spectrograph: on 
October 18, 1994 and on October 29, 1998. The third one oc- 
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Fig. 7. Flare observed on February 9, 2007 in Budapest. The bi- 
nary and spot model is removed from the observed light curves. 
Light curves from up to down are made using B, V, Rc, Ic filters. 
Dashed lines show zero level for each passband (an arbitrary 
shift was applied for better visibility). Arrows show the times of 
the post-flare events. 



curred on February 9, 2007, and was detected in BV(RI)c pho- 
tometric bandpasses. 



3.3.1. Spectroscopic events 

The 1994 flare was visible in Ha, HB and He I D3 lines. During 
the flare, the Ha equivalent widths increased (see Fig. 0), in- 
dicating stronger activity in the chromosphere. On the dynamic 
spectrum (Fig. [8] left panel) as well as on the individual spectra 
(FigHJ left panel) enhancement in the blue wing of the Ha line is 
observed since the beginning of the observations. The flare erup- 
tion started at phase 0.51, reached its maximum at phase 0.58 
and lasted for more than 2.5 hours, at least until the end of the 
observations. At the time of the flare HB is also increased, and 
asymmetry is seen in the red wing at all phases (Fig. [9] middle 
panel). As the indicator of the explosive event He I D3 appears 
as prominent emission line (Fig. [9] right panel) at ph ases 0.54- 
0.58. A very similar event is reported on LO Peg by Eibe et al. 
(Il999h . 

The weaker flare from October 29, 1998 shows increased 
Ha emission, and the HB level is also higher, but only a marginal 
He I D3 emission is seen (see Fig.QT]in the online version). This 
weaker flare was not accompanied by a (large) prominence, the 
Ha line is quite symmetric, except before and during the flare 
at phases 0.06-0.14 (Fig. [8] middle panel). The Ha line outside 
the small flare is the same as in the quiescent Ha state (Fig. [8] 
right panel). The HB line is again increased in the red wing at all 
phases. 

The increase of the blue wing of the Ha line is probably 
the signature of prominence formation, i.e. emerging material in 
the magnetic loop with an upflow of ~40km s (see Fig. [TOlin 
the online version) which forms a prominence. After the flare 
the Ha line is again symmetric, since most of the material is 
heated up and leave the loops during the fast phase of the flare 
eruption. The red asymmetry of the HB line could originate from 
the downfiow of cooling material along the flux loops. 
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Fig. 8. Dynamic Ha spectra from October 18 and 23, 1994, October 29 and 28, 1998 (from left to right). Crosses on the right show 
the time of observations, circles on the plot and stars mark the radial velocity curve of the primary and secondary, respectively. The 
first two images on the left show two flares, the third one is from a quiescent state of the binary. All plots show phases between 
0.1-0.9. The left plot shows observations from October 18, 1994 except the last two points from October 23, 1994 after phase 0.75. 
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Fig. 9. Ha, HjS and Na D lines around the time of the flare event in October 18, 1994 between phases 0.37-0.74. The flare occurred 
between phases -0.5-0.7 (see scale on the right). Note that the scale of the H/3 plot is different because of the higher noise in the blue 
part of the spectrum. Diamonds on the Ha spectra show the place of the secondary component calculated from the radial velocity 



The Ha line of the secondary star is much better visible on 
the spectra made at the quiescent state of the binary (Fig. fT2lin 
the online version), which is partly due to the better signal-to- 
noise ratio of these spectra, but also the less active primary at 
this time does not wash out the weaker signal of the secondary. 
During the observations of the quiescent phase, in the beginning, 
some red asymmetry is seen again in the HB line, indicating that 
chromospheric flows are present on the star in most of the time. 

3.3.2. Photometric flare 

The light curve of the flare on February 2007 is plotted in Fig. [7] 
The system at the peak of the flare is 0.6 magnitudes brighter in B 
than the quiescent light, and the flare is well seen also in Iq filter, 



indicating a very strong outburst. The event could be observed 
for 3.72 hours, which is about one-third of the rotation. Using 
photometric data, unfortunately no details can be determined for 
the position of the flare. It is possible that the flare lasted longer 
than ~4 hours, but the flaring region moved out of view because 
of the rotation. Since both components are active, the flare could 
occur in either of the two components. 

Assuming that the whole light curve of the flare was ob- 
served, we can make an estimation of the flare energies. The en- 
ergy emitted by the flare in B, V,Rc and Ic bands are 1.43e+35, 
7.58e+34, 1.34e+35 and 1.31e+35 ergs, assuming that the flare 
occurred on the primary. 

If we suppose that the flare occurred on the much fainter 
secondary, we need to know its magnitudes for calculating flare 
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energy. According to CI97, the secondary component is at least 
3 magnitudes fainter than the primary. If we assume solar metal- 
licity of the system, and a mass o f 0.2 1M Q from our results fo r 
the secondary, using the tables of |VandenBerg & ClerrJ (|2003), 
we get B = 16.1, V = 14.4, R c = 13.1, I c = 11.5 magnitudes as 
an upper limit of the secondary component's brightness. Using 
these values and taking into account the uncertainties in stellar 
temperatures and radii (see Table [2]i the resulting flare energies 
satisfactorily agree with the values above, that we determined 
supposing the flare occured on the primary. For th e details of 
flare energy calculating, see e.g. lKovari et alJd2007|) . 

After the peak of the flare, four smaller, post-flare eruptions 
were observed. These oscillations could be the result of recon- 
nection between the moving footpoints of the expanding mag- 
netic loop causing the flare and the interm ittent magnetic field 
(see iKovari et alJl2007l: lAttrill et all 120071) or individual addi- 
tional flares. However, the quasip eriodicity of the events makes 
the first scenario more likely (cf. lMitra-Kraev et al.ll2005l) . 



4. Discussion 

Our measurements show that V405 And is a very active RS CVn 
binary. Strong flare events were observed both in photometry 
and spectroscopy (see Section [3~3l l. Photometric data shows that 
the stellar surface is covered by large spotted regions, which are 
quite stable during a timescale of the observations (~650 day), 
i.e., the maximum spot visibility (spotted light curve minima) is 
near the same phase. The eruption on September 30, 2007 hap- 
pened at this position, and probably the one on September 29, 
2007 near to that. The large flare on October 18, 1994 also ap- 
pears at the same phase. These evidences point towards the pres- 
ence of an active nest that has a very long lifetime, more than 10 
years. The phase of the eruptive events is near to the secondary 
minimum, so the feature is probably locked to the binary orbit. 
Such l ong-lived active regions hav e been observed e.g. on WY 
Cnc by lKozhevnikova et al] (120071) . 

According to CI97, log(L x /L bo i) for V405 And i s -3.1+0.14 
which fits well the rotation act ivity relations of iJames et"ai] 
( 2000) and Pizz olato et ail d2003l) 



We compared the masses and radii of V405 And (Table [2]) 
with the existing measured values of dwarf stars with masses 
below 0.8 M B in Fi g. [T31 The theoretical curve for 5 Gyr from 
iBaraffe et all £l998) is also plotted using r eff , = 5780 and 
A^boi.Q = 4.72. Note that at this part of the theoretical mass-radius 
diagram there is not much difference between isochrones be- 
longing to different ages; low mass stars evolve slowly. It is well 
seen that stars less massive than ~O.32M , which is suggested 
as limit of full convection, fit well the isochrones. Among these 
low-mass stars are the secondary components of the binary from 
NGC 1647 and of V405And. However, with increasing mass 
the correlation gets looser. In Fig. [13] there are three systems 
in which the masses of the components are very differen t: T- 
Lyrl-17236, the binary from NGC 1647 dHebb et al.ll2006h and 
V405 And. T-Lyrl-17236 rotates rather slowly with a period of 
P m t - 8.4c/, whereas the other two systems are very fast rotators 
(V405 And: P mt 0.465c/, star from NGC 1647:P rot * 0.619c/). 
The radii of the primary components of the two fast rotating 
systems are well above the theoretically predicted radius val- 
ues, especially V405 And is much larger than expected. On the 
other hand, both components of the slow rotator T-Lyrl-17236 
fit rather well the isochron e. 

iMullan & MacDonaldl (120011) found that the radii of low- 
mass M dwarfs, concerning stellar structure codes, are larger 
than expected, and that could be due to their strong magnetic 
fields, which m ay push the m ass limit of full convection towards 
lower masses. Riba| d2006l) compared theory with observation 
on masses and radii of low-mass stars (both binaries and single 
stars). From the high precision results of double lined eclips- 
ing binaries he found that theory predicts smaller radii by about 
10% (or higher temperatures by a bout 5%) for stars in th e mass 
range between 0.4-0. 8M . Finallv. lChabrier et al.l(l2007l) earned 
out evolutionary calculations on low-mass stars taking into ac- 
count magnetic activity, and demonstrated that magnetic fields 
alter the evolution, and that the spot coverage yields larger radii 
and smaller effective temperature. However, the radius of the pri- 
mary co mponent of V405 And is even larger comparing to the 
results of Cha brier et al.l (120071) . Future modellings with mag- 
netic fields are necessary to map the area of M dwarfs in the 
mass-radius plane. An interesting question is, what would be the 
effect of the magnetic interaction between a very active, fully 
convective star with its close, also active component in a binary 
system, so that the radius of the more massive component be- 
comes even larger than the recent theory predicts, while at the 
same time the radius of the fully convective star agrees with the 
prediction. 

At present only two systems are known whose primaries are 
about 2.5 times more massive than the secondaries, and that one 
component is above while the other is below the limit of the full 
convection. The primary of V405 And deviates most from the 
expected radius of its mass from the known sample. On the other 
hand, the secondary of V405 And would be fully convective even 
in the prese nce of its strong magnet i c field according to the cal- 
culations of IMullan & M acDonaldl d200ll Fig. 1). It would be 
a challenge to model evolution of such binaries whose compo- 
nents are on the two sides of the mass limit of full convection. 
Of the two known examples, V405 And is the one with its V ~ 
1 1 mag. brightness which is well observable showing all the sig- 
natures of strong magnetic activity like spots and flares. We rec- 
ommend this binary to the attention of future investigators. 
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Fig. 10. Dynamic Ha spectrum with the flare from October 18, 1994 showing the excess emission compared to the average of all 
spectra. The spectra are shifted with their actual radial velocity. 
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Fig. 11. He, H/3 and Na D lines around the time of the flare event in October 29, 1998 between phases 0.56-0.30. The flare occurred 
between phases ~0. 05-0. 15 (see scale on the right). Scales and notation is the same as in Fig. [9] 
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Fig. 12. Hot, H/3 and Na D lines from the quiescent state of V405 And in October 28, 1998. The contribution of the secondary can 
be well seen around the Ha line. Scales and notation is the same as in Fig. [9] 



